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Abstract 

We consider the electroweak radiative corrections in the MSSM and study the purely Yukawa 
contributions 0(a L yf p y^ q ), p + q = L, where y^,6 are the top and bottom quark Yukawa couplings. 
We show that these corrections can be computed in a gaugeless limit of the MSSM where they are 
under Renormalization Group control. As an application, we present explicit results for various 
ILC and LHC processes valid at all orders in the loop expansion and at leading order in the large 
logarithms that arise at high energy. 

PACS numbers: 11.30.Pb, 12.15.Lk, 11.15.Bt 



I. INTRODUCTION 



i] 



The Minimal Supersymmetric Standard Model (MSSM) []J is one of the major theoret- 
ical laboratories where physics beyond the Standard Model can be tested at quantitative 

n 

level. Radiative corrections in the MSSM are under control at all orders [2j, at least in 
processes where a perturbative treatment is adequate. Of course, the large number of free 
parameters of the MSSM remains a main practical problem. The comparison between the- 
ory and experiment is difficult to perform in a model independent way and with no specific 
assumptions. 

Some simplification in the treatment of radiative corrections could be quite useful to 
disentangle individual and different physical mechanisms. High energy expansions of virtual 
effects are a natural candidate in such direction 3]. Indeed, both at the Large Hadron 
Collider (LHC) and at the International Linear Collider (ILC) the experimental setup will 
be characterized by a large available energy compared to the typical mass scales, at least in 
favorable (but not at all exotic) scenarios. 

High energy expansions are based on the simple observation that asymptotic scattering 
amplitudes are dominated by large radiative corrections growing with energy. For brevity, we 
shall indicate collectively all these contributions as Sudakov logarithms (SL), although only 
a part of them deserves this terminology. A detailed quantitative analysis of the Sudakov 
effects is by now very well understood both in the Standard Model 0] and in the MSSM 5|. 

At one loop and in 2 — * 2 processes, SL can be classified in four very different classes §|. 
We concentrate on the electroweak sector and discuss the four categories emphasizing in 



particular the status of their theoretical control. The detailed application of this formalism 



to realistic physical problems can be found in |7[ for ILC processes and in 



a 



for LHC 



processes. 

The first class is that of Universal SL which are infrared contributions. The and 
Z° gauge boson masses act as infrared regulators in the loop diagrams and we can observe 
enhanced corrections at high energy with leading form at one loop ~ a log 2 (s/M^), where 
\fs is the center of mass energy and a = e 2 /4ir. These mass singularities are called universal 
because they receive a contribution from each external particle (initial or final) independently 
of the diffusion details. Also, the weight of the squared logarithm is a simple combination 
of the gauge group quantum number of the external particles. At L loops, the leading 
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contributions of this kind are ~ a L \og 2L (s/M^ v ) with subleading corrections involving all 
powers smaller than 2L. The universal SL are the most relevant from the numerical point 
of view. For this reason there have been large efforts to compute them at higher orders or 
(working at fixed order in a) to extend the calculation to additional subleading terms. The 
theoretical control of the universal SL is coded in the hard evolution equations. Detailed 
applications at the two loop level in the Standard Model can be found in Q]. Explicit 
diagrammatic calculations have also been accomplished 

In the MSSM, the existing calculations are at the one loop level, with the exception of a 
certain number of resummation proposals valid at all orders in the perturbative expansion 
and at NLO in the logarithms [13j . 

The next class we consider is that of Renormalization Group (RG) SL. These are com- 
plementary to the universal ones. Indeed, RG-SL are effects with an ultraviolet origin and 
arise as a consequence of the running of the coupling constants. Being a short distance effect 
they are expected to be independent on the long distance physics. At L loops they have the 
leading behavior ~ a L \og L (s/M^). Since there is only one logarithm per loop and not two 
as in the universal case, the RG Sudakov contributions are NLO at one loop, NNLO at two 
loop and so on. Thus, although they are completely under control from the theoretical point 
of view, they are not expected to be the dominant contribution in any reasonable limit. 

The third class is composed of the Angular SL. At one loop, they are contributions coming 
from Standard Model box diagrams and with asymptotic behavior ~ alog(s/M^) log \x/s\ 
where x = t or u, the Mandelstam variables involving the scattering angle. The Stan- 
dard Model contributions have been investigated at two loops [11]. They do not receive 
supersymmetric corrections and we shall regard them as known without further discussion. 

The fourth and last class is the most interesting for our purposes and consists of the 
Yukawa SL. These contributions are present only when heavy quarks or their SUSY partners 
are produced in the final state (or belong to the initial state, at LHC). More precisely, 
Yukawa SL are effects with characteristic factors m^/M^y where m t ^ are the top and bottom 
quark masses Q|. These factors are of course Yukawa couplings in disguise. Again, at 
leading order, we have one logarithm for each loop. In the Standard Model, Yukawa SL are 
numerically non negligible, but not so interesting for the same reasons we mentioned about 
RG-SL. However, in the MSSM the Yukawa couplings can be enhanced in certain regions of 
the parameter space and Yukawa SL permit interesting phenomenological analyses. These 



3 



regions are characterized by large values of the mixing parameter tan/3. This important 
parameter is singled out in the Yukawa SL weights and it is possible to fix bounds on it from 
the Sudakov expansion of the observables of several processes [14 1. 

The theoretical control over Yukawa SL is not very satisfactory as it stands. In the 
Standard Model, they have been often neglected by working in processes involving only light 
fermions ^(|. In the MSSM, where the top quark physics programme is very promising, the 
existing calculations are limited to the one loop level. 

On the other hand, the one loop analysis suggests that a more detailed investigation 



should be possible because of some remarkable features [1J, [JJ|. First, the Yukawa SL 
turn out to be correction factors associated to the external particles and not to the specific 
details of each process. This is similar to what happened with the universal SL. Second, 
there is only one logarithm per loop and this seems to be the marker of a ultraviolet effect. 
A posteriori, it is strange that there are RG-SL related to the gauge coupling running, but 
nothing analogous in the Yukawa sector where additional independent couplings are present. 

In this paper, we fill these gaps and explain the above features. We show that Yukawa 
SL are indeed a short distance effect which is governed by Renormalization Group equations 
in the gaugeless limit of the MSSM, i.e. in the Yukawa sector. Contrary to the RG-SL, we 
have motivations to pursue a calculation of Yukawa SL at higher order. We anticipated the 
reason and we repeat it here to emphasize its importance. There are regions in the MSSM 
parameter space where the Yukawa couplings can be enhanced. If we are in such a region, 
the Yukawa SL can give large corrections at one loop, comparable to universal SL, although 
they are subleading in a formal logarithmic expansion. For the purposes of a high precision 
measurement, the determination of the next perturbative SL contributions acquire therefore 
a substantial relevance. 



II. YUKAWA ENHANCED ELECTROWEAK CORRECTIONS IN THE MSSM 
AT HIGH ENERGY 

The detailed structure of Yukawa enhanced Sudakov corrections in the MSSM is illus- 
trated by two sample processes shown in Fig. The first example, e + e~ — > ti, is a typical 
s-channel neutral process and describes top quark pair production at ILC. The second ex- 
ample is instead a charged process, ud — > tb, and is one of the three partonic processes 



4 



FIG. 1: Two sample processes relevant for ILC and LHC respectively. In both processes the leading 
Yukawa SL are confined in the dashed circles representing the final dressed interaction vertex. 



which permit single top quark production at LHC [12|. In the Feynman diagrams we have 
drawn the final vertex by means of a dashed circle standing for a particular set of radiative 
corrections. These are precisely the electroweak virtual exchanges that give rise to purely 
Yukawa SL, i.e. all contributions that at the L loop level and at high energy correct the 
tree level amplitude by factors that take the form 

where m t ^ are the top and bottom masses, s is the squared cm. energy, and M is a 
typical process mass scale. We emphasize that there are also SL with factors like above, 
but with p + q < L. We shall not discuss these contributions which mix the infrared and 
— structure of the MSSM. lu part.eu.ar, the papers Q resuru SL at NLO order 
(that is including terms ~ log 2L_1 s) and determine at L > 1 the contribution of the form 
Eq. (P) with p + q = 1 and r = 2L — 1. Our calculation extends their result by independent 
additional terms. We deal with the corrections having p + q = L which we can prove to 
have a clean ultraviolet origin. As a consequence, the exponent r turns out to be bounded 
as r < L in this case. Our study will also be further limited by the condition r = L, i.e. we 
shall calculate the leading contributions in the logarithmic expansion. 

Since the combined power of factors m t ^ is maximal in Eq. it is impossible to generate 
Yukawa SL in other 1PI parts. For instance, gauge boson self energies will always have too 
many gauge couplings. Box diagrams will have suppressed Yukawa couplings to initial state 
light fermions. As in the above examples, similar diagrams occur in other processes like in 
particular those related to Fig. (JTJ) by supersymmetry, as discussed in jis| . 

Due to the counting of couplings, the Yukawa effects in Eq. (JTJ) can be computed in an 
essentially gaugeless limit of the MSSM. In particular, the radiative corrections shown in 



Fig. ([1} require the calculation of a vertex diagram with two on-shell chiral and antichiral 
fields and an external classical vector field, neutral or charged according to the process. 

In this gaugeless limit, the only relevant piece of the MSSM lagrangian is the Yukawa 
superpotential 

W = y t t (W° u - bH+) + y b b (tH d - bH° d ). (2) 

The chiral fields t, b contain the right handed antitop and antibottom, the SU(2) doublet 
Q = (t, b) is composed of the chiral fields t and b containing the left handed top and bottom, 
the SU(2) doublets H u = (H+,H®) and Hd = (H^H^) contain the various Higgs and 
Higgsino fields of the MSSM. The couplings y t} b have the tree level value (g is the SU(2) 
gauge coupling and Mw is the W boson mass) 

g m t I g m b 1 

Vt ~ v^M^sin/3' Vb ^/2M w cos/3' U 

in terms of the conventional vacuum alignment angle j3. 

About soft breaking terms, we shall see in a moment that the corrections Eq. arise 
as a short distance effect and are thus quite independent on the soft breaking Lagrangian. 

After these preliminary discussion and definitions, we lift the discussion from the examples 
of Fig. (Q) to the specific universal source of Yukawa enhanced logarithms. This is the basic 




vertex diagram V — > shown in Fig. (J2J), where V is a charged or neutral gauge 

superfield, depending on the specific final state Let us denote by g the relevant gauge 

coupling. We want to compute the corrections 0(g yl P yl q h L ), p + q = L, to the one particle 
irreducible (1PI) vertex in Fig. (J2J). We denote the subset of such corrections as t[J\Q) 
where Q is a large momentum entering on the gauge field line and where the two matter fields 
are on-shell. Notice that we shall keep working in manifest supersymmetric formalism by 
employing superfields. Indeed, perturbation theory will be trivial at the leading logarithmic 
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accuracy with all non trivial effects being encoded in the Renormalization Group (RG) 
equations. 

The fact that rg°(Q) can be computed in a gaugeless limit with only external classical 
gauge fields and not internal quantum ones is very important and has deep consequences. 
Indeed, under these conditions, the large momentum behavior of Y^\Q) is governed by 



Renormalization Group (RG) evolution equations 



hard evolution equations jl(| . Technically, in the Callan-Symanzik equation for T\ v - 1 (Q) the 
mass insertion term is irrelevant because there are no infrared singularities. As a conse- 
quence, one is left with the much simpler Renormalization Group equation which we can 
solve. This means that the large logarithms of Yukawa origin are actually a genuine short 
distance effect that can be addressed in the deep Euclidean region. In particular, they are 
expected to be independent at leading order on the various mass terms appearing in £ so ft- 
The RG evolution equation for is 

d , a i.. .. a d , a ,.. .. \ 9 .. (.. .. \\ t^(V) ( Q 
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171 ] and does not require the so-called 



<y)i 



where (3 t ,b are the Yukawa /3-functions, 7^ = 7, + jj, and ji is the anomalous dimension of 
$j. All RG functions are computed in the gaugeless limit g — > 0. The scale \i is the RG 
subtraction mass. In deriving Eq. (j3J) we neglect systematically all non leading terms in the 
y t) b expansion like gauge boson self energy corrections. 

In the next Section, we shall analyze the consequences of Eq. (j3J) at the leading logarithmic 
accuracy. 



III. RENORMALIZATION GROUP ANALYSIS 



As is well known, the general solution of the evolution equation Eq. (jlj) is 



r 



(V) 



-,yt,Vb 
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Vt 
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n 



Vb 



■ Vb 
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X 



(5) 



where Ut,b{Q/^) are the running Yukawa couplings according to the MSSM Yukawa (3- 
functions. Eq. (JSJ) expresses the renormalized vertex at the scale Q in terms of its value 
at the initial scale /i. In principle, we should need the perturbative expansion of rj^(l, y). 
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However, if we are only interested in the leading logarithms 0(K log we can argue from 
Eq. (jSJ) that they appear in a single coefficient multiplying the tree level value of the vertex 

where depends only on the one loop terms in the anomalous dimensions and /3-functions. 

To give an explicit expression for Cy it is convenient to define the coefficients *, 7^ b 
appearing at one loop in the anomalous dimensions (a tt b = Vtb/i^ 71 )) 



(1) t Olt (1) bOib 

7 * = 7 * 4^ + 7 * + 



as well as the following sums over the external i, j fields 



(7) 



t (i) t . (i) t 

>y. . — ry: ' 4- 'V. 
/ 1 j It ' 3 ' 



b (1) 6 , (1) 6 



(8) 



After some algebra, the coefficient Cj,- turns out to be given by the following compact formula 







a b (/i 2 ) 






[MQ 2 )\ 



with < 



^(67^-7?,), 
(67?/ -ihl 



70 
1 



(9) 
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where the running couplings a t ,b(Q 2 ) are consistently evaluated at leading logarithmic ac- 
curacy (i.e. solving the one loop RG equations). 

The practical evaluation of Eq. (JHJ) requires (a) the various anomalous dimensions at 
one loop, (b) the solution of the one-loop RG evolution equations for the running Yukawa 
couplings. We discuss separately each issue and also provide a perturbative expansion of 
Eq. whose validity we shall discuss in the applications. 



A. Renormalization Group functions 

The RG functions (anomalous dimensions and /^-functions) can be found in the litera- 
ture [l8[ for general N = 1 gauge models. In the gaugeless limit they reduce at the one 
loop level to the following simple results that we present according to the field numbering 
($1, . . . , $ 8 ) = ti, H®, bi-, H+, 1>r, Hj, H®). The first three anomalous dimensions are 

7i = 2(f) 2 , i2 = ^y 2 {y 2 t +yll ^ = ^yvl (io) 
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The other ones are 74 = 72, 75 = 73. Also, 76 and 77 are obtained from 71 and 73 with the 
replacement y t <-> yb- Finally 73 = 77- 

The Yukawa coupling /3-functions for the two couplings y t and yb are 

A = |p K^o + + $ ' A = (^2 K^cf + 3)y b 2 + yt] • (11) 

B. Running a t}b (Q 2 ) 

From the one loop /3-functions we determine the running Yukawa couplings. In terms of 
s = Q 2 and setting N c = 3 the evolution equations read 

^ ^(Qa t + a b ), (12) 



d log s An 
da b a b 
d log s 4tt 

with given initial conditions at the scale fi. These equations can be solved analytically in 
implicit form but the exact solution is not particularly enlightening. We provide full details 
in Appendix [X] Beside the exact solution, it is interesting to have also its perturbative 
expansion permitting to discuss the convergence of the loop expansion. At second order, we 
have 

where in the r.h.s. we denote a. t .b = a tj b{^ 2 )- The result for db{Q 2 ) / atbi.^ 2 ) is obtained by 
exchanging a t <-> «6- The expansion Eq. (fT3j) can be substituted in Eq. (jUJ). Expanding 
again in powers of the logarithms we obtain the leading logarithm approximation for the 
Sudakov correction Cjj (no sum over repeated indices) 



% = 1 - lQ g ^(^«* + 7l» + (14) 



+ W ^ [«? 7&(7& - 12) + of 7&(T& " 12) + 

Eqs. (JSJ) and its two loop expansion Eqs. (J 14)) are the main result of this paper. The one 

re first line) is in perfect agreement with the results already published 
and coming from explicit diagram calculations in component fields 
followed by high energy expansions. In the next Section, we shall discuss in details the 
consequences of our results by considering specific applications. 



loop approximation ft 
in the literature 



IV. APPLICATIONS AND DISCUSSION 



A. Preliminary perturbative analysis 

We begin with a discussion of the perturbative result at two loop order Eq. (|14j) . This 
will be useful in comparing the one loop approximation with higher order corrections. We 
denote L tJb = ^ log%. From Eq. (JHD, we can obtain the following list of specific cases 



tR t* R : 


: c u = 


1 - 


2L t 


~ ^\ - 


LtLb, 


tL t* L 


: c 22 = 


1 - 


L t - 


-L b - 5 - 
2 


(L 2 t + Ll 


b L b* L 


: c 44 = 


C22, 








tL b* L 


: c 24 = 


C22, 








bR b* R 


: c 55 = 


1 - 


2L b 


- ~ 


■ LtLb, 


K h- 


: c 55 = 


1 - 




2 * 


~ — L+Lh, 
2 




: c 77 = 


1 - 


3L b 


V- 
2 Lb 


- —LfLh. 
2 



The meaning of Eq. is as follows. Let us consider for instance the first line which 
reports the expression of Cu. This is the correction factor for the vertex which has as a final 
state the chiral field (and its conjugate) whose fermionic component is the right handed top 
quark. By supersymmetry, the same correction factor is also obtained if we consider in the 
final state various combinations of the scalar partners, for instance tnt* R . In such a diagonal 
case, the vector V is necessarily a neutral one, 7 or Z. Also, by supersymmetry, from the 
real superfield V we can also take its gaugino component. 

Notice that the physical charged Higgs boson is a mixture of the charged components of 
H u and H^. A proper two loop result for the vertex with this field would require a one loop 
treatment of the mixing that we defer to later work. The vertices with final quark superfields 
(describing production of quark pairs, squark pairs or quark-squark combinations) are correct 
as they stand, since mixing is irrelevant when it appears only in the internal lines. 

The above logarithmic expansions Eqs. (fTK|) can be resummed in closed form at least in 
the single Yukawa coupling limit. For instance, if we are in a point of MSSM parameters 
where we can make the approximation y b y t , then it is possible to repeat the derivation 
leading to and exploit the exact solution for the running a t (Q 2 ). The final result is 
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remarkably simple 

y b = : Clj = (1 - 6L t yh/ 12 . (16) 

As an example, in the case of final states t R i,R or t^L, we find (in agreement with the L t 
terms in Eq. (fTHj) ) 

t R t* R : en = (1 - 6L t ) 1/3 = 1 - 2L t - 4L 2 + • • • , (17) 
*L ^ : c 22 = (1 - 6L t ) 1/6 = l-L t -h 2 t +--- . 

B. Discussion of the full result 

Now, we turn to consider the numerical relevance of the computed effects according to 
Eq. (jnj) with both Yukawa couplings being active. We concentrate on the same processes 
shown in Fig. (0), i.e. 

ILC e+e-^ f a J a , f = t,b 

LHC u d -> t L b L . 

In both cases the polarization of the initial state affects only the tree level amplitude and is 
factored in the correction. The chirality index a determine the polarization of the final state 
in the ILC process. For brevity, we consider only a final state with asymptotic vanishing 
helicity. The final state of the LHC process has a = L because of W boson exchange. In 
both cases, the correction factor in the cross section is simply (cy) 2 . 

The explicit values of the exponents rf^ required for the evaluation of Eq. (0) are sum- 
marized in the following table where, for clarity, we replace the pair with the fermion 
in the associated chiral field. As we remarked, by supersymmetry the same corrections are 
obtained with the sfermion components. 



*, 3 


rf 


if 


tR tR 


12/35 


-2/35 


t L t L 


1/7 


1/7 


b L b L 


1/7 


1/7 


b R b R 


-2/35 


12/35 



The initial condition for a t (Q 2 ), o>b{Q 2 ) is fixed as follows. We choose a scale /i and define 
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at that scale 



( 



yt_ 

An 



An 



) 



) 



2 



2 




(1 +COt 2 /3 eff ) 



(1 + tan 2 /3 e fr), 



(19) 



where the parameter tan /3 e ff is a scale dependent effective mixing angle. At tree level, this 
is the conventional mixing angle. Beyond tree level, it is just a convenient parametrization 
of the initial values of the top and bottom Yukawa couplings. For simplicity, we shall denote 
in the following discussion (3 = (3 e g. 

Our results are shown in the three Figs. (013) for the ILC processes. The left hand side 
of each figure shows the (Yukawa) Sudakov correction to the cross section at one loop and 
two loop level. We also show a line with the label exact standing for the evaluation of the all 
order Eq. (|9~|) . The initial condition is obtained from Eq. ()19j) with tan/3 = 40. The curves 
are shown as functions of \fs with the RG scale is fixed at the arbitrary value fi = 100 GeV, 
somewhat between and m t . As usual, the choice of the scale /i cannot be fixed at the 
leading order in the logarithmic expansion. 

In all the considered cases the two loop correction is practically equivalent to the exact 
result. Instead, the correction to the one loop approximation is appreciable. In the right 
hand side of the three figures, we show this difference as a function of \fs for two reference 
values of tan (3 = 2, 40. 

The value of the extra correction is in all cases at the level of 1-2 %, which can be visible 
at ILC. The only exception is the process e + e~ — > bn bn at small tan f3 due to the very small 
value of yb in that case. 

About the considered LHC process, the equality C24 = C22 implies that the correction 
factor is identical to that occurring in e + e _ — ► t L t L . However, there are important differences 
deserving a couple of comments. First, the available energy for the partonic process will be 
in general smaller than at ILC reducing the overall effect. Second, it will be non trivial to 
identify it with the invariant mass of the final state tJ>L due to various systematic corrections 
(for a realistic analysis of these effects in tt production, see the second reference in 0). 
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V. SUMMARY AND CONCLUSIONS 



This paper has been devoted to the analysis of a specific subset of the radiative corrections 
affecting physical processes in the MSSM. In the framework of a high energy expansion we 
have concentrated on large logarithmic corrections to the invariant amplitudes A that take 
the following asymptotic form at L loops 

A = A^ll + ^t^ log L 4)> ( 20 ) 

V L>1 ^ ) 



p=0 



We have shown how to compute the coefficients {d^p} for various specific processes at ILC 
or LHC by solving the Renormalization Group equation governing the relevant 1PI parts 
from which the above correction originates. The final recipe is quite simple and universal. 
It explains various features of the above correction already noted in the literature as a 
byproduct of explicit component calculations. 

In the Standard Model, the coefficients {d^p} are fixed numerical constants. Instead, 
in the MSSM, they can vary by large amounts as the MSSM parameter space is explored. 
In benchmark scenarios where one of the two relevant Yukawa couplings is large, we can 
obtain sizable corrections with relevant two loop contributions. The next order corrections 
are rather small in all the presented examples. 

We remark that such two loop contributions are difficult to obtain in a straightforward 
approach that would consist in evaluating the relevant Feynman diagrams at generic kine- 
matical points, followed by analytic high energy expansions. For instance the two loop 
contribution is NNLO in the logarithmic expansion whereas the very powerful tools pre- 
sented recently in work only at NLO as they stand. Instead, the calculation of {d^ p } is 
quite easy within the framework of the Renormalization Group. 

The main point of our paper has been precisely to emphasize that the Yukawa Sudakov 
corrections can be actually computed in a gaugeless limit where they arise not as mass 
singularities, but as ultraviolet effects. This observation relies on old results about the 
asymptotic properties of quantum field theory form factors and combines them with the 
modern phenomenological interest toward large Yukawa coupling regions of the vast MSSM 

rn 

parameter space. A typical example of such scenario is the conventional SPS4 point [20]. Of 
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course, this kind of strategy is well known for static quantities where large logarithms also 



appear. A nice example is the recent calculation (21 1 of gluino lifetime and branching ratios 
in Split Supersymmetric models where the large corrections involve \og(rh/rrig), where m is 
the squark and slepton scale and m-g is the gluino mass. In that case, the large logarithms 
controlled by the top quark Yukawa coupling are also resummed separately. 

The aim of our analysis has been that of showing that similar simple RG treatments 
are possible in the context of high energy expansions and in particular for the Yukawa SL 
whose existence is somewhat independent of the more involved genuinely Sudakov double 
logarithms that require more sophisticated tools. From a numerical point of view, our 
analysis has shown that the higher order SL effect could be of the percent size at ILC, 
in which case it would become visible. Less optimistic conclusions would apply to the 
considered LHC processes for which the effect appears to be definitely small. Note, however, 
that this fact represents in any case a valuable information, telling us that in this sector the 
available one loop expansions are certainly sufficient. 

Various extensions of the present analysis are possible. For instance, two physically 
interesting examples at LHC are it pair production via the partonic processes gg — *> it, 
qq — ► it and the single top quark production mechanisms that we have not discussed, 
i.e. bu — > td and bg — > tW . They admit similar corrections localized in specific three 
point vertices 0, |^ to which the present formalism can be applied to obtain higher order 
contributions. 
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APPENDIX A: SOLUTION OF THE COUPLED RG EQUATIONS 



The equations to be solved are 

a t {t) = a t {Qa t + a b ), (Al) 
a b (t) = a b (6a b + a t ), (A2) 

where we have introduced the variable t = 4- log(s/so) that we shall call time in the following. 
This dynamical system admits the family of integral curves 

(a t - a b ) 7 = C (a t a b ) 6 , (A3) 



where C is a constant that can be computed at the initial point t = t . Exploiting Eq. (|A3j) . 
we can write a tyb parametrically in terms of the ratio R = a t /ct b . It is necessary to treat 
separately the two cases R(to) > 1 and R(to) < 1. If R(to) = 1, then C = and R = 1 at all 
times. In this case we have simply the single coupling RG solution at = a b = at(0) /(l — 7t). 
It is enough to work out the case R(to) > 1, since the other is obtained by the swap a t <-> a b . 
The parametric equations are 

a t = C- 1/5 R- 1/5 (R-l) 7/5 } (A4) 
a b = C-^R-^tR-l) 7 '^ 

where the ratio R satisfies 

R(t) = 5C- 1/5 R~ 1/5 (R - 1) 12/5 . (A5) 

We immediately see that R increases as t increases. In practice there is an exploding time 
at which R — ► oo. This happens at unphysical large energies for realistic initial values. 
Solving the separable differential equation Eq. (jA5|) . we find the (implicit) solution for the 
ratio R = R(t) in terms of its initial value Rq = R(to) 

F(R)-F(R o ) = C- 1 ' 5 (t-t ), (A6) 

where the monotonically increasing function F(R) is 

'R^ S (R + 1) l fn ^3/5 _ /3 4 8 



(i?- l) 7 /5 



(A7) 
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Thus, an accurate and quite simple recipe to solve the initial system is the following three 
steps algorithm: (i) compute C from the initial values of a^&, (ii) at each time find numeri- 
cally the unique zero of Eq. (jA6|) . (iii) replace the root R(t) in Eqs. (|A4|) . This semi-analytic 
procedure is superior to any discrete step integration scheme, like for instance Runge-Kutta 
methods. Indeed, from Eq. (jA6|) we derive immediately exact properties of the solution. An 
example is the (unphysical) exploding time which is predicted by Eq. (|A6|) to be 

/ t + r^(F f(r )) f - 1 r(i/5)r(8/5) 
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